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Abstract: A selfcontained proof of the KAM theorem in the Thhring model is discussed. 
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I I shall particularize the Eliasson method, [E], for KAM tori to a special model, of great interest, whose 

■ relevance for the KAM problem was pointed out by Thirring, [T] (see [G] for a short discussion of the model) . 
T-H I The idea of exposing Eliasson's method through simple particular cases appears in [V], where results of the 

pH . type of the ones discussed here, and more general ones, are announced. 

^ I The connection between the methods of [E] and the tree expansions in the renornialization group approaches 

■ to quantum field theory and many body theory can be found also in [G[. The connection between the tree 
0^ I expansions and the breakdown of invariant tori is discussed in [PV]. 

The Thirring model is a system of rotators interacting via a potential. It is described by the haniiltonian 

" ^ I (see [G] for a motivation of the name): 

Cn; lj-^A-A + ef{d) (1) 

O . 2 

I where J is the (diagonal) matrix of the inertia moments, A — {Ai, . . . , Ai) G are their angular momenta 

■ and d = (ai, ... ,0;;) G T' are the angles describing their positions: the matrix J will be supposed non 
fT^ I singular; but we only suppose that miuj^i^. .^; Jj = Jg > 0, and no assumption is made on the size of 
0^ ' the twist rate T — min J^"^: the results will be uniform in T (hence the name "twistless": this is not a 

\ contradiction with the necessity of a twist rate in the general problems, see problems 1, 16, 17 in §5.11 of 

[G2], and [G]). We suppose / to be an even trigonometric polynomial of degree TV: 

'V ■ 

O; /(«) = X! fucosiy-d, fp = f-iy (2) 

, 0<\P\<N 
^ ■ 

• • , We shall consider a "rotation vector" cug = (wi, . . ■ ,uJi) G i?' verifying a strong diophantine property (see, 

^ ' however, the final comments) with dophantine constants Co,t, 7,0; this means that: 

H' 1) coMo-^^i > i^?r", Q^u^z' 

■ " " ' 2) min ICoMo ■v\-Y\> 7"+^ if n < 0, < \v\ < (7"+^)"^ ^ 

0>p>n 

and it is easy to see that the strongly diophantine vectors have full measure in i?' if 7 > 1 and c are fixed 
and if r is fixed t > I — 1: we take 7 = 2, c = 3 for simplicity; note that 2) is empty if n > — 3 or p < n + 3. 
We shall set Aq — Jujq. And a special example can be the model fo{d) — JoWQ(cosai + cos(q:i + 0:2)). 

We look for an e -analytic family of motions starting at d — and having the form: 

A^ Ao + H{ujot;e), d = Q^t + h{ujQt; e) (4) 

with H{ijj ; e), h{i/j ; e) analytic in ip G T' and in e close to 0. We shall prove that such functions exist and 
are analytic for \ limjjj \ < ^ for \e\ < Eq with: 

eoi=6Jo-ic2/o^^+'e^^e«^ (5) 
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where b, c are I -dependent positive constants, fo = max^; \fff\. This means that the set A = AQ + H{tp ; e), a = 
tp + h{'ip ; e) deseribed as ip varies in is, for e small enough, an invariant torus for ( 1 ), whieh is run quasi 
periodically with angular velocity vector wq • It is a family of invariant tori coinciding, for e = 0, with the 
torus A = Ao, a = tp G T'. One recognizes a version of the KAM theorem. The proof that follows simplifies 
the one reported in [G]. 

Supposing Jo = Ji < J2 the uniformity in J2 (i.e. what we call the twistless property) implies that the same 
£0 can be used as an estimate of the radius of convergence in e of the power series describing the KAM tori 
with rotation vector ljq — (^1,^2) in the system {2Ji^Al + LU2A2 + e/o(ai,a2), which is one of the most 
studied hamiltonian systems. The estimate can be improved. Note that a careful analysis of the proof of 

the KAM theorem also shows the uniformity in the twist rate in the case (1). 

Calling H'-'^\tp),h'-'^\'ip ) the fc-th order coefficients of the Taylor expansion of H,h in powers of £ and 
writing the equation of motion as a = + J~^A and A = —sdsf{d) we get immediately recursion relations 
for Namely ojq -dhf^ = Jj^Hf^ and: 

Q,.dHf = - p^^^^«.•C-^„-77('I5o^)•E*^^^^'^°^('^o^) (6) 

|75i|>0 

where the ^* denotes summation over the integers kj>l with: X]i=i X]j=i kj = k — 1. 
The trigonometric polynomial /i'*^-' {tp ) will be completely determined (if possible at all) by requiring it to 
have average over %p , (note that i?*^*^) has to have zero average over ip ). For k = 1 one easily finds: 

h^'Ki^) = li^^e--^" (7) 

Suppose that lv-^\tp) is a trigonometric polynomial of degree < kN , odd in t, for 1 < fc < A'q. Then we 
see immediately that the r.h.s. of (6) is odd in t. This means that the r.h.s. of (6) has zero average in 
t, hence in ip , and the second of (6) can be solved for k = ko. It yields an even function H^'^^^ip) which 
is defined up to a constant which, however, must be taken such that H^''°\tp) has zero average, to make 
Co ■dhj'^^ = JJ^H^^^ soluble. Hence the equation for /i^*^^ can be solved (because the r.h.s. has zero average) 

and its solution is a trigonometric polynomial in ip , odd if h^''^ is determined by imposing that its average 
over 'ip vanishes. 

Hence the (7) provide an algorithm to evaluate a formal power series solution to our problem. It has been 
remarked, [E],[V], see also [G], that (6) yields a diagrammatic expansion of h^''\ We simply "iterate" it until 
only h^^\ given by (7), appears. 

Let 1? be a tree diagram: it will consist of a family of "lines" {i.e. segments) numbered from 1 to A; arranged 
to form a (rooted) tree diagram as in the figure: 

jl 

vi ve 



Vu 

fig. 1: A tree diagram 1? with m„j, =2,mi,j =2,m„2 =3,m^3 =2,my^ = 2 and m = 12, ]^ m„! = 2"* .6, and some decorations. 
The line numbers, distinguishing the lines, are not shown. 



root 
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To each vertex v we attach a "mode label" iyy G Z, \uy \ < N and to each branch leading to v we attach a 
"branch label" jy = 1, . . . ,1. The order of the diagram will be fc = number of vertices = number of branches 
(the tree root will not be regarded as a vertex). 

We imagine that all the diagram lines have the same length (even though they are drawn with arbitrary 
length in fig. 1). A group acts on the set of diagrams, generated by the permutations of the subdiagrams 
having the same vertex as root. Two diagrams that can be superposed by the action of a transformation 
of the group will be regarded as identical (recall however that the diagram lines are numbered, i.e. are 
regarded as distinct, and the supcrpositon has to be such that all the decorations of the diagram match). 
Trees diagrams are regarded as partially ordered sets of vertices (or lines) with a minimal element given by 
the root (or the root line). We shall imagine that each branch carries also an arrow pointing to the root 
("gravity" direction, opposite to the order). 

We define the "momentum" entering v as i'{v) = X)^>^, I'w If from a vertex v emerge mi lines carrying a 
label j = 1, m2 lines carrying j = 2, . . ., it follows that (6) can be rewritten: 

m _ iv-* -TT (iJ~v^)3„./g„nLi(H)r° f . 

with the sum nmning over the diagrams t9 of order k and with ^{vq) = u] and the combinatorics can 
be checked from (6), by taking into account that we regard the diagram lines as all different (to fix the 
factorials). 

The * recalls that the diagram d can and will be supposed such that ^{v) ^ for all u G (by the above 
remarked parity properties). Note that (8) is implied by the corresponding (6.23) of [G]: one can check that 
the two formulae coincide (by summing over what in [G] are called the "fruit values"). The theory in [G] is 

in fact a little more general, although it is really applied to the same Thirring model. 

There are other diagrams, however, which we would like to eliminate. They are the diagrams with nodes 
v' ,v, with v' < V, not necessarily nearest neighbours, such that u{v) = v{v'). Such diagrams are, according 

to Eliasson's terminology, in resonance. 

To see this, for the purpose of illustration, we shall first restrict the sum in (8) to a sum over diagrams such 
that: 

I) v{v) ^ Q if u is any vertex. 

II) i/{v) ^ i^W) for all pairs of comparable vertices v' (not necessarily next to each other in the diagram 
order, however), with v >Vo. 

There are at most 2'^''k\ diagrams and the labels j can be at most l'^ while the Vy can be chosen in a number 
of ways bounded by {2N + 1)"= < (3^^)'*=. Therefore, if /o = max;; \fv\- 

\hi^h < Ar(3iV)"=22'=Z'=Ml^Ar2feniaxTT(CoWo • j^{v))~^ < {foC^Jo^)''N^^+^^''+\2l3^)''M (9) 
"^0 ^ 

where the maximum is over the diagrams 'd verifying the I), II) above. Hence the whole problem is reduced 
to estimating the maximum with M. 

Let C{q) = (SN'^q^)^^ , where r is the diophantine constant in (3): then Col^o • i^l > 5({q) if < < q: 
we can say that e is "g-singular" if \uo ■ y\ < ({q). Then the following (extension) of a lemma by 
Brijuno, see [P], holds for diagrams of degree k verifying I), II) above. Fixed g > 1 let N{k) be the number 
of q-resonant harmonics, among the vertex momenta h'{v) of the vertices v G i). Then N{k) < 2-|. 

Assuming the above claim true (we shall not use it outside the present heuristic argument), we fix an 
exponentially decreasing sequence 7"; the choice 7 = 2 recommends itself. The number of g = 7~"-singular 
harmonics which are not q = 7~^"'"'"^^-singular is bounded by 2fc7~", (being trivially bounded by the number 
of 7 ""-singular harmonics!). Hence: 

^ —00 

IlTCoSrW- n C(7^"-^^)-^'=^"<exp4fc^2-"log(3.2^("-)iV^)=M (10) 

V u u \ ;j n<-l n>l 

therefore the series for the approximation to ft, V , that we are considering because of the extra restriction in 
the sum (8), has radius of convergence in e bounded below. 
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However there are resonant diagrams. The key remark is that they cancel almost exactly. The reason is 

very simple: imagine to detach from a diagram the subdiagram with first node v. Then attach it to all 
the remaining vertices w > v',w G ■d/'d^- We obtain a family of diagrams whose contributions to h'^^^ differ 
because: 

1) some of the branches above v' have changed total momentum by the amount ?(«): this means that some 
of the denominators (w • i'(w))~^ have become (w • ^{v) ± e)^"^ if £ = cDq • v{v)] and: 

2) because there is one of the vertex factors which changes by taking successively the values v^j, j being the 
line label of the line leading to v, and w G'&/'d2 is the vertex to which such line is reattached. 

Hence if a; ■ = e = we would build in this resummation a quantity proportional to: ^ — v{v) — y{v') 
which is zero, because v{v) = v{v') means that the sum of the v^'s vanishes. Since u ■ v = e ^ Q we can 
expect to see a sum of order e^, if we sum as well on a overall change of sign of the values (which sum 
up to 0). 

But this can be true only if £ <C cD • i?', for any line momentum V of a line in i?/t?2- If the latter property 
is not true this means that Q ■ v' \s small and that there are many vertices in '&/'&2 of order of the amount 
needed to create a momentum with small divisors of order e. 

Examining carefully the proof of Brjiuno's lemma one sees that such extreme case would be esscintially also 
treatable. Therefore the problem is to show that the two regimes just envisaged (and their "combinations") 
do exhaust all possibilities. 

Such problems are very common in renormalization theory and arc called "overlapping divergences" . Their 
systematic analysis is made through the renormalization group methods. We argue here that EUiasson's 
method can be interpreted in the same way. 

The above introduced diagrams will play the role of Feynman diagrams; and they will be plagued by 
overlapping divergences. They will therefore be collected into another family of graphs, that we shall call 
trees, on which the bounds are easy. The {Q ■ are the propagators, in our analogy. 

We fix an scaling parameter 7, which we take 7 = 2 for consistency with (3), and we also define tS = ColSq: 
it is an adimensional frequency. Then we say that a propagator (w • i')"^ jg scale n if 2"~^ < jw • i7| < 2", 
for n < 0, and we set n = 1 if 1 < |d; ■ 

Proceeding as in quantum field theory, see [G3] , given a diagram ?9 we can attach a scale label to each line 
v'v in (8) (with v' being the vertex preceding v): it is equal to n if n is the scale of the line propagator. Note 
that the labels thus attached to a diagram are uniquely determined by the diagram: they will have only the 

function of helping to visualize the orders of magnitude of the various diagram lines. 

Looking at such labels we identify the connected clusters T of vertices that are linked by a continuous path 
of lines with the same scale label ut or a higher one. We shall say that the cluster T has scale tit- 

Among the clusters we consider the ones with the property that there is only one diagram line entering 
them and only one exiting and both carry the same momentum. Here we use that the diagram lines carry 
an arrow pointing to the root: this gives a meaning to the words "incoming" and "outgoing" . 

If V is one such cluster we denote Ay the incoming line: the line scale n = n\^, is smaller than the smallest 
scale n' = ny of the lines inside V. We call wi the vertex into which the line Ay ends, inside V. We say 
that such aV is a resonance if the number of lines contained inVis< E2~"^, where n = n\y, and E,e 
are defined by: E = 2~^^N~'^, e = r^^. We shall say that n\y is the resonance scale. 

Let us consider a diagram ^ and its clusters. We wish to estimate the number Nn of lines with scale n < 
in it, assuming Nn > 0. 

Denoting T a cluster of scale n let m-p be the number of resonances of scale n contained in T {i.e. with 
incoming lines of scale n), we have the following inequality, valid for any diagram ^?: 

Nn<^^+ E (-1 + ^t) (11) 
T, nT=n 

with E = N~^2~^^, e = t~^. This is a version of Brjuno's lemma: a proof is in appendix. 
Consider a diagram 1)^ we define the family T{'d^) generated by i!)^ as follows. Given a resonance V of 
we detach the part of above Ay and attach it successively to the points w £ V, where V is the set of 
vertices of V (including the endpoint wi of Ay contained in V) outside the resonances contained in V. Note 
that all the lines A in have a scale > ny. 

For each resonance V of i?^ we shall call My the number of vertices in V. To the just defined set of diagrams 
we add the diagrams obtained by reversing simoultaneously the signs of the vertex modes i^j, for w €V: the 
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change of sign is performed independently for the various resonant clusters. This defines a family of Y\ "^My 
diagrams that we call The mimbcr Y[ 2Afy will be boimdcd by cxp^^ 2Mv < e^*^. 

It is important to note that the definition of resonance is such that the above operation (of shift of the 
vertex to which the line entering V is attached) does not change too much the scales of the diagram lines 
inside the resonances: the reason is simply that inside a resonance of scale n the number of lines is not very 
large being < iV„ =£2-"^". 

Let A be a line, in a cluster T, contained inside the resonances V = Vi C V2 C . . . of scales n = ni > n2 > ■ ■ 
then the shifting of the lines Ay; can cause at most a change in the size of the propagator of A by at most 
2"i + + . . . < 2"+^ 

Since the number of lines inside V is smaller than Nn the quantity tJ-ux of A has the form Lj-i^ + a\uj ■ v\y 



if i7? is the momentum of the line A "inside the resonance V. 



it is the sum of all the vertex modes of 



the vertices preceding A in the sense of the line arrows, but contained in V; and (j\ = 0, ±1. 
Therefore not only |w • > 2"+^ (because is a sum of < Nn vertex modes, so that < NNn) but 
cJ ■ is "in the middle" of the diadic interval containing it and by (3) docs not get out of it if we add a 
quantity bounded by 2"^^ (like a\Oj ■ vxy). Hence no line changes scale as "d varies in J^{'d'^), if cDq verifies 
(3). 

This implies, by the strong diophantine hypothesis on loq, (3), that the resonant clusters of the diagrams in 
) all contain the same sets of lines, and the same lines go in or out of each resonance ( although they 
are attached to generally distinct vertices inside the resonances: the identity of the lines is here defined by 
the number label that each of them carries in -d^). Furthermore the resonance scales and the scales of the 
resonant clusters, and of all the lines, do not change. 

Let 1?^ be a diagram not in J^{'&^) and construct J^{'&'^), etc. We define a collection {^(i?')}i=i,2,... of 
pairwisc disjoint families of diagrams. We shall sum all the contributions to /i^*^^ coming from the individual 
members of each family. This is the Eliasson's resummation. 

We call Ev the quantity w ■ v\y associated with the resonance V. If A is a line with both extremes in V 

wc can imagine to write the quantity lu ■ i/x as cJ • + (TxSv, with ax = 0, ±1. Since |a; • > 2"^~^ we 
see that the product of the propagators is holomorphic in sy for \sv\ < 2"^~^. ^ While £y varies in such 
complex disk the quantity |w • does not become smaller than 2"^~^ — 2 2"^~^ > 2"^~^. Note the main 
point here: the quantity 2"''^'' will usually be ^ 2"^v~^ which is the value ey actually can reach in every 
diagram in J^['&^); this can be exploited in applying the maximum priciple, as done below. 
It follows that, calling nx the scale of the line A in t?^, each of the Yl 2Mv < e^'^ products of propagators 
of the members of the family !F{d^) can be bounded above by Ylx 2~^("^~^) = 2^"^ Ylx 2~^"^, if regarded as 
a function of the quantities ey = w • i^Xv^ for \sv\ < 2"^~^, associated with the resonant clusters V. This 
even holds if the sy are regarded as independent complex parameters. 

By construction it is clear that the sum of the Y[ 2Mv < e^*^ terms, giving the contribution to h!^^^ from the 
trees in J^id^), vanishes to second order in the sy parameters (by the approximate cancellation discussed 
above). Hence by the maximum principle, and recalling that each of the scalar products in (8) can be 
bounded by N'^, we can bound the contribution from the family by: 



fc! V Jo ; 



n<0 



■2nJV„ 



n n 

n<OT, nT=n i=l 



2(ra-ni+3) 



(12) 



where: 

1) Nn is the number of propagators of scale n m. 'd^ (n = 1 does not appear as • i7| > 1 in such cases), 

2) the first square bracket is the bound on the product of individual elements in the family T{'9^) times the 
bound e?'^ on their number, 

3) The second term is the part coming from the maximum principle, applied to bound the resummations, 
and is explained as follows. 

i) the dependence on the variables ey = Si relative to resonances Vi dT with scale uxy = n is holomorphic 
for for \ei\ < 2"*~'^ if m = nvi, provided rij > n + 3 (see above). 

ii) the resummation says that the dependence on the Si's has a second order zero in each. Hence the maximum 
principle tells us that we can improve the bound given by the first factor in (12) by the product of factors 



' In fact > 2"+^ because F is a resonance; therefore Iw-j^aI > 2"+^ — 2"+-^ > 2"+^ so that ny > n + 3. 

On the other hand note that |a; • > 2"^-^ - 2"+^ so that \uJ ■ + axSyl > 2"^-^ - 2"+^ - 2"^-^ > 
2"v-i _ 2"v-3 > 2»v-2^ fQj. jg^i < 2"-''-3. 
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{\si \ 2~"<+^)2 if m > n + 3. If Tii < n + 3 we cannot gain anything: but since the contribution to the bound 
from such terms in (12) is > 1 wc can leave them in it to simplify the notation, (of course this means that 
the gain factor can be important only when <C 1). 

Hence substituting (11) into (12) we see that the tut is taken away by the first factor in 2^"2~^"*, while the 
remaining 2~^"* arc compensated by the —1 before the +mT in (11), taken from the factors with T = Vi, 
(note that there are always enough — I's). 

Hence the product (12) is bounded by: 

iiV(Co'^o"VoiV2)'^e2'^24'=26'=[] 2-6"'=^"2" < (^3^ 

n 

with: Bo = 2iOe2C^/oiV2exp[iV(22+2^"' log 2) E^iP2-P""']. 

To sum over the trees we note that, fixed i? the collection of clusters is fixed. Therefore we only have to 
multiply (13) by the number of diagram shapes for (< 2^*^^!), by the number of ways of attaching mode 
labels, (< (3iV)"=), so that we can bound by (5). 

Comments. 

The strong diophantine condition is quite unpleasant as it seems to put an extra requirement on luq: I 
think that in fact such condition is not necessary. Given loq verifying the first of (3) with some constant 
Cq and some r we can imagine to define Co = this leaves the first of (3) still valid. One should then 

note that the < < (2"+^)""^ ^ will imply that if < < (2"+^)"'^ ^ the numbers \iJ ■ v\ are spaced 
by at least (2(2'^2"+3)-'^"')-'^ = 2"+^. Hence we can find a sequence 7„ such that 1 < 7„2-" < 2 and 
||w • z7| — 7p| > 2"+^ for > J9 > n. Defining a propagator x"^ to have "scale n" if 7„_i < |x| < 7„ it should 
be possible to perform the proof without the second assumption in (3). This will be interesting to check if 
it really works as it would be fully constructive, given verifying the first of (3), and as it would cover the 
most studied cases like the case / = 2 and Q = (u;i,u;2) = (l)'^) with uj a quadratic irrational, or a number 
with uniformly bounded continued fraction entries. It has also the conceptual advantage that the sequence 
of scales 7„ is not "prescribed a priori" but it is determined by the arithmetic properties of the rotation 
vector ujQ. 



Appendix Al: Resonant Siegel-Brjuno bound. 

Calling the number of non resonant lines carrying a scale label < n. We shall prove first that A^* < 

2k{E2-^-^)-^ liN* > 0. 

If has the root line with scale > n then calling 192, ■ • ■ , '&m the subdiagrams of emerging from the 
first vertex of d and with % > £;2-^" fines, it is N*{d) = N*{'di) + . . . + N*{'drn) and the statement is 
inductively implied from its validity for k' < k provided it is true that N*{'&) = if A: < E2~^^, which is is 
certainly the case if E is chosen as in (11).^ 

In the other case it is N* < 1 + Y^^i and if m = the statement is trivial, or if m > 2 the 

statement is again inductively implied by its validity for k' < k. 

If m = 1 we once more have a trivial case unless the order ki of t^i is fci > — ^E2~^^. Finally, and this 
is the real problem as the analysis of a few examples shows, we claim that in the latter case the root line is 
either a resonance or it has scale > n. 

Accepting the last statement it will be: N*{d) = I + N* {di) = 1 + N* {•&{) + ... + N*{-d'^,), with -d'j being 
the m' subdiagrams emerging from the first node of iS'i with orders kj > E2^^'^: this is so because the root 
line of will not contribute its unit to N*{'di). Going once more through the analysis the only non trivial 
case is if m' = 1 and in that case N{'&[) = N*{'d^'i) + . . . + N{-&"m"), etc, until we reach a trivial case or a 
diagram of order < fc — i_E2~"^. 

It remains to check that if fci > A; — ^E2~"^ then the root line of i?i has scale > n, unless it is entering a 
resonance. 

Note that if A; < _B2-"^ it is, for all momenta u of the lines, < NE2-"^, i.e. \uj ■ u\ > {NE2-"^)-'^ = 2» 2" so that there 
are no clusters T with nx = n and N* = 0. The choice E = N~^2~^^ is convenient: but this, as well as the whole lemma, 
remains true if 3 is replaced by any number larger than 1. The choice of 3 is made only to simplify some of the arguments 
based on the resonance concept. 
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Suppose that the root Une of i?i is not entering a resonance. Note that \ij ■ t'(wo)| < 2", \t3 ■ < 2", 

if vo,vi arc the first vertices of i? and respectively. Hence S = \(oj ■ (i7(wo) — iy(vi))\ < 2 2" and the 
diophantine assumption imphes that |i'(t'o) — i^ivi)] > (2 2")""^ , or i^{vo) = i^{vi). The latter case being 
discarded as fc — fci < 2~"^ (and we are not considering the resonances), it follows that k — ki< \E 2~"^ 
is inconsistent: it would in fact imply that v{vq) — v{vi) is a sum of /c — fci vertex modes and therefore 
\v{vq) — i^{v\)\ < ^NE2~"^ hence (5 > 2^ 2" which is contradictory with the above opposite inequality. 
A similar, far easier, induction can be used to prove that if N* > then the number p of clusters of scale 
n verifies the bound p < k {E2~'^")~^ — 1. Thus (11) is proved. 

Remark: the above argument is a minor adaptation of Brjiuno's proof of Siegel's theorem, as remarkably 
exposed by Poschel, [P] . 

Acknowledgements: I am indebted to J. Bricmont and M. Vittot for many clarifying discussions; and in 
particular to S. Miracle-Sole for the same reasons and for his interest and suggestions.. 
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